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abstract 

We present a light curve analysis and radial velocity study of KOI-74, an eclipsing A star + 
white dwarf binary with a 5.2 day orbit. Aside from new spectroscopy covering the orbit of 
the system, we used 212 days of publicly available Kepler observations and present the first 
complete light curve fitting to these data, modelling the eclipses and transits, ellipsoidal mod- 
ulation, reflection, and Doppler beaming. Markov Chain Monte Carlo simulations are used to 
determine the system parameters and uncertainty estimates. Our results are in agreement with 
earlier studies, except that we find an inclination of 87.0 ± 0.4°, which is significantly lower 
than the previously published value. The altered inclination leads to different values for the 
relative radii of the two stars and therefore also the mass ratio deduced from the ellipsoidal 
modulations seen in this system. We find that the mass ratio derived from the radial veloc- 
ity amplitude = 0. 104 ± 0.004) disagrees with that derived from the ellipsoidal modulation 
= 0.052 ± 0.004 assuming corotation). This was found before, but with our smaller inclina- 
tion, the discrepancy is even larger than previously reported. Accounting for the rapid rotation 
of the A-star, instead of assuming corotation with the binary orbit, is found to increase the 
discrepancy even further by lowering the mass ratio to ^ = 0.047 ± 0.004. These results indi- 
cate that one has to be extremely careful in using the amplitude of an ellipsoidal modulation 
signal in a close binary to determine the mass ratio, when a proof of corotation is not firmly 
established. The same problem could arise whenever an ellipsoidal modulation amplitude is 
used to derive the mass of a planet orbiting a host star that is not in corotation with the planet's 
orbit. 

The radial velocities that can be inferred from the detected Doppler beaming in the light 
curve are found to be in agreement with our spectroscopic radial velocity determination. We 
also report the first measurement of R0mer delay in a light curve of a compact binary. This 
delay amounts to — 56±17s and is consistent with the mass ratio derived from the radial 
velocity amplitude. The firm establishment of this mass ratio at ^ = 0. 104 ± 0.004 leaves little 
doubt that the companion of KOI-74 is a low mass white dwarf. 

Key words: binaries: close - binaries: eclipsing - stars: individual (KOI-74). 



* Partly based on observations made with the Mercator Telescope, oper- 
ated by the Flemish Community, and the William Herschel Telescope, both 



located at the Spanish Observatorio del Roque de los Muchachos of the 

Instituto de Astrofisica de Canarias. 

t E-mail: steven.bloemen@ster.kuleuven.be 



©2011 RAS 



2 S. Bloemen et al. 



1 INTRODUCTION 

The primary science goal of the Kepler Mission is the detection of 
Earth-Uke exoplanets, but its highly accurate photometric observa- 
tions also reveal hundreds of eclipsing binary stars (^Slawson et al. 
l201lh and are well suited for t he study of stellar varia bility at un- 
precedentedly low amplitudes JPebosscher et al.ll201lh . In this pa- 
per, we present an analysis of 212 d of Kepler data of the eclipsing 
binary KOI-74 (KIC 6889235), and spectra taken at different or- 
bital phases. 

The system consists of a main- sequence A star primary and 
a less massive companion. Its light curve shows deeper eclipses 
than transit^ which implies that the companion is hotter than the 
primary. There is a clear asymmetric ellipsoidal modulation pat- 
tern in which the flux maximu m after the t ransit s is larger than 
the maximum after the eclipses. TRowc et aP bOld) suggested that 
the asymmetry in the ellipsoidal modulation be due to a star spot. 
Van Kerkwijk et al. (2010) instead attributed it to Doppler beam- 
ing. This effect is caused by the stars' radial velocities which shift 
the spectrum, modulate the photon emission rate and beam the pho- 
tons somewhat in the direction of motion. The effect was , as far as 
we are aware, first discussed in'S hakura & Postnovl (Il987h and first 
observed by Maxted et al. (2000). Its expect ed detection in Kepler 
lig ht curves was s ugges ted and discussed bv lLoeb & Gaudil (l2003h 
and lZucker et al.l (l2007h . The detection of Doppler beaming in Ke- 
pler light curves has led to t he discovery of seve ral non-eclipsing 
short-period binary systems Jpaider et al.l |201 ih and it has been 
sho\ yn that it can also be observed in planeta r y systems (see 
e.^. | Mazeh & Faidej |201QI : IShporer et all 1201 ll) . iBloemen et aP 



ID detected Doppler beaming in the Kepler light curve of 
KPD 1946+4340 and presented the first comparison between a ra- 
dial velocity amplitude derived from Doppler beaming with a spec- 
troscopic value. In the case of KPD 1946-1-4340, the results were 
found to be consistent. For KOI-74, spe ctroscopic radial velocit y 
measurements were recently presented bv lEhrenreich et aP (l201lh . 
In this paper, we present independent spectroscopic radial veloc- 
ity measurements which we compare with the photometric radial 
velocity amplitude prediction. 

Earlier analyses of the K epler light curve of KOI - 74 ar e 
presented in iRowe et al.l bold) and Ivan Kerkwiik et al ] 12013). 
ikowe et aP (|201Q|) measured the mass ratio of the system from the 
amplitude of the ellipsoidal modulation and found a companion 
mass of 0.02-0.11 M0. Van Kerkwijk et al. (2010) b uilt on these 
result s (e.g. they used the inclination value report by iRowe et aP 
[20i3) but claimed that the mass ratio of the system could not be 
determined reliably from the ellipsoidal modulation amplitude. In- 
stead, they used the radial velocity information from the Doppler 
beaming signal to derive a companion mass of 0.22 ± 0.03 M©. 
They concluded that the companion has to be a low mass white 
dwarf and showed that the system properties are in good agreement 
with a binary that has undergone a phase of stable Roche lobe over- 
flow from the more massive star to the less massive star. This puts 
KOI-74 in an evolutionary stage th at follows on that o f systems 
such as WASP J0247-2515, which iMaxted et all (1201 ih recently 
identified as a binary consisting of an A- star and a red giant core 
stripped from its envelope. Up to now, 4 close binaries consisting of 
a white dwarf and a main sequence star of spectral type A or F have 



^ We use the terms 'transit' and 'eclipse' to indicate, respectively, the oc- 
cultation of the A- star by the compact object and the occupation of the 
compact object by the A-star half an orbit later. 
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Figure 1. Fit to the spectral lines of KOI-74 using the average HERMES 
spectrum (after shifting the spectra t o the rest frame of th e primary) and the 
solar metaUicity synthetic spectra of lMunari et al.l ( l2005h . The uncertainties 
on the parameters indicated on the figure only reflect the formal errors on 
the fit. 



been found in Kepler data ("Rowe et al. I|2010l:l van Kerkwiik et al.l 
2010; Carter et al. 2011; Breton et al. 20lll). 

We remodel the Kepler light curve of KOI-74, adding an addi- 
tional 175 d of data compared to the previous studies and perform 
Markov Chain Monte Carlo (MCMC) simulations to explore the 
uncertainty on the derived system parameters. We will revisit the 
issue of the true mass ratio of the system, using input from the 
light curve analysis (Doppler beaming, ellipsoidal modulation and 
R0mer delay) and spectroscopy. 



2 SPECTROSCOPY 

To determine the spectral type, the rotational velocity and the ra- 
dial velocity amplitude of the primary, we obtained 46 high resolu- 
tion {R ~ 85000) echelle spectra using the HERMES spectrograph 
(Raskin etal. 2011.) at the 1.2-m Mercator telescope (La Palma, 
Canary Islands). These were reduced using the standard instru- 
ment specific data reduction pipeline. Additionally, we obtained 29 
spectra with the ISIS spectrograph mounted on the 4.2-m William 
Herschel Telescope (La Palma, Canary Islands). We averaged the 
HERMES spectra, after shifting out the primary's radial velocity. 
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Figure 2. Radial velocity curve of the primary of KOI-74, measured by 
fitting a Gaussian profile to the mean profile of three Balmer lines (see text 
for details). The radial velocity measurements from ISIS@WHT spectra are 
represented by black circles, the measurements from HERMES @Mercator 
spectra by green squares, folded on the orbital period. The systemic velocity 
derived from both datasets has been subtracted. The error bars are scaled to 
deliver a unit %^ per degree of freedom. 



From the average spectrum, we measured the rotational velocity 
of the prima ry vsin/ = 1 64±9kms~^, using the Fourier method 
presented in lGravl (Il992h , from the Mg II XA4SI line. Compared 
to other techniques such as fitting broadened synthetic spectra, the 
Fourier method has the advantage of bei ng rather unsensitive to 
other line broadening mechanisms (see e.g. lSimon-Diaz & Herrerol 
l2007h . The Mg II M481 line is a doublet with components at 
4481.126 A an d 4481.325 A, wh ich leads to an overestimation of 
vsin/ (see e.g. lRoveretal.|[2002h . By applying the same method 
to a synthetic spectrum with a comparable rotational broadening, 
we find that we can expect our vsin/ measurement to be over- 
estimated by about 10km s~^ due to the double nature of the 
line. Given this overestimation, our re sult is in line with t he val - 
ues mentioned in t he Note in proof of Ivan Kerkwiik et al ] (l201Qh . 



150kms"^ and in lEhrenreich et all tOllk 145 zb 5kms~^ If the 
primary were in corotation with the binary orbit, we would expect 
vsin / ~ 25kms~^ The rotational velocity measured from spec- 
troscopy thus implies that the primary is not in corotation but is 
instead a fast rotator. We have adopted v sin / = 150±10kms-i for 
the analysis presented in this paper. We also performed a spectral 
fit to the Balmer li nes and Mg II X44S 1 using the solar metallicity 
synthetic spectra of iMunari et al.l (l2005h and assuming a rotational 
velocity of 150km s~^. The fits to the spectral lines are shown in 
Fig. [T] The small emission feature in the core of the Ha and Hp 
lines does not originate from any of t he two binary compone nts 
but is caused by sky emission (see also lEhrenreich et aDl201ll) . It 
is absent in the ISIS spectra, for which we could perform a back- 
ground subtraction during the data reduction. We find that the pri- 
mary has a gravity of logg ~ 4.27 and an effective temperature 
of Tgff ~ 9500K. These results are in agreeme nt with the s pectra l 
type AlV and T^ff = 9400± 150K derived bv lRowe etaP ilOld) . 
from which they inferred a primary mass of Mi = 2.2±O.2M0. 
We adopted T^ff = 9 500 ± 250 K for our analysis. 

We have used two techniques to derive the radial velocity 
amplitude of the A- star. The only metal line that is clearly de- 
tected in both the HERMES and ISIS spectra, is the Mg II line 
at 4481 A. We measured the radial velocities from that line by fit- 
ting a Gaussian profile to it. We also measured the radial veloci- 
ties from the Balmer lines at 4102A, 4340A and 486lA simulta- 
neously, by first normalising the sp ectra using low order splines 
(as described in 'Pap ics et aP 120121) . and then fitting a Gaussian 
to the core of the mean profile obtained by Least Squares Decon- 



Table 1. Radial velocity amplitudes (Ki)md systemic velocities (y) of KOI- 
74, measured from ISIS@WHT and HERMES @Mercator spectra. The un- 
certainties that are given are scaled to get a unit per degree of freedom. 
Initial reduced is given in the last column. Two techniques have been 
used: a Gaussian fit to the Mg II line at 4481 Aand a Gaussian fit to the 
mean profile of the Balmer lines at 4102A, 4340A and 486lA(see text for 
details). 



Instrument 


Line(s) 


Ki (kms-^) 


y(kms ^) 


5Cinit,reduced 


ISIS 


Balmer 


15.8±0.4 


-50.9 ±0.3 


1.2 




Mg II 


15.4±0.7 


-49.3 ±0.5 


1.6 


HERMES 


Balmer 


14.9 ±0.4 


-47.1 ±0.3 


0.5 




Mg II 


16.5±1.3 


-51.0±1.0 


2.3 


Weighted mean (adopted) 


15.4±0.3 


-49.1 ±0.2 





volution (lDonatietal.ll 19971) . The measured radial velocity ampli- 
tudes and systemic velocities, assuming a circular orbit, are given 
in Table [T] The uncertainties on individual data points have been 
scaled to deliver a unit per degree of freedom. The initial x^, 
before rescaling the uncertainties, is given in the table. The radial 
velocities measured from the Balmer lines are shown in Fig. |2l 
fo lded on the orbital period using the Kepler ephemeris as given 
in lRowe et al.l (|201Q|) . Some spectra were taken in bad seeing con- 
ditions, which is reflected by the large error bars on a few of the 
data points. We have adopted the weighted mean velocity ampli- 
tude of our measurements, Ki = 15.4 ± 0.3km s~^, for the analysis 
presented in this paper. The weighted mean systemic velocity, is 
y = 49.1 ± 0.2km s~^. Our radial velocity amplitudes are in agree- 
ment with (but slightly lower than) Ehrenreich et al. (2011)'s re- 
sult of 18.2 ± 1.7kms~^. We provide the radial velocities we have 
measured from our spectra in electronic form with this paper. 



3 KEPLER PHOTOMETRY 

The Kepler data from QO (quarter 0), Ql, Q2 and Q3 were retrieved 
from the public archivq^- The data span 229 d, resulting in a dataset 
of 212 d of observations excluding the gaps. The data of QO, Ql 
and the first two months of Q2 are long cadence data (30 m integra- 
tions), the last 27 d of Q2 and the entire Q3 dataset (86 d) are taken 
in short cadence mode (1 m integrations). The light curve (see Fig. 
[3] in this paper for a phase folded version) shows cl ear eclipses, 
transits, el lipsoidal modulation and D oppler beaming. iRowe et aP 
( 201 3) and lvan Kerkwiik et al.l ( l201(ib already presented models for 
the light curve. We remodelled the light curve with the L CURVE 
code written by TRM, using mor e data, and find a lower orbital 
inclination than lRowe et"aD (l201Qh (/ = 88.8 ± 0.5°). This finding 
also has consequences for som e of the system parameters derived 
by van Kerkwiik et al] (l201Qh . as they adopted the inclination of 
Irowc et al.l (12010 ) for their analysis. 

Below we discuss the various steps of our binary light curve 
modelling. 



3.1 Detrending of the Kepler data 

We extracted the Kepler data from the pixel data, rather than using 
the pipeline output. We first normalised the data by dividing out an 



^ Publicly released Kepler data can be downloaded from 
http :// archive .stsci . edu/kepler/ ^ 
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Figure 3. The Kepler QO, Ql, Q2 and Q3 data of KOI-74, folded on the orbital period and binned into 2000 phase bins. The black curve shows a typical model 
fit to the data. The residuals of the data points after subtracting the model are shown in the bottom panel, binned into 2000 phase bins (green) and 200 phase 
bins (black). 



initial model fit, to be able to determine the trends in the data. Next, 
we detrended the data by dividing out exponentials fitted to the two 
obvious instrumental decays during Q2 and Q3. We removed 101 
long cadence points from the Q2 data, as well as 1299 short cadence 
datapoints from the Q3 data, which were too severely affected by 
instrumental effects to allow for accurate detrending. We then di- 
vided out second order polynomials fitted to each segment of the 
light curve. While fitting the polynomials, we flagged 184 points as 
> 4-a outliers, in an iterative fashion, leaving 166 836 datapoints. 
We then folded the initial orbital fit back in. Finally, we binned the 
short cadence data out of eclipse to long cadence data since there 
are no binary signatures to be expected outside eclipse that require 
the time precision offered by short cadence. This way, the number 
of datapoints is reduced to 28 086. 



Initial MCMC runs using this dataset (see Section 13.41 ) 
showed, as expected, that there is no correlation between the 
adopted ephemerides {T{) and P) and the other free parameters. We 
therefore decided to create a phasefolded version of the four months 
of short cadance data to speed up the simulations. We rebinned the 
light curve into 1800 s bins out of eclipse and 30 s bins in eclipse, 
to be left with 1717 data points. 



We assume that there is no background contamination in the 
Kepler data. We provide our final datasets in electronic form with 
this paper. 



3.2 What can we learn from the Kepler light curve? 

Fitting a light curve of a binary that shows a variety of effects (such 
as eclipses, ellipsoidal modulation and Doppler beaming) allows 
one to gain a lot of information on the binary's parameters. We first 
explore the potential before we discuss the actual light curve fits in 
Section[3H 

A first useful constraint comes from the measured fractional 
transit depth, which reveals the ratio of the radii of the two stars: 

AFtransit/F = {R2/Rlf- (1) 

Given R2/R1 from the transit depth and Ti from the spectral fit 
(Section|2]), the eclipse depth provides information about T2. Using 
blackbody approximations and neglecting the contribution of the 
companion out of eclipse, one finds: 

F ^ Fi Qxp{hv/kT2)-\ \RiJ ' 

The eclipse duration te gives as a function of the inclina- 
tion / (assuming spherical stars): 

(Rl/a)^ = sin^ {%te/P) sin^/ + cos^/ (3) 

teu sselll [ml) with P the orbital period and a the separation be- 
tween the two binary components. 
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In case of corotation of the star with the binary orbit, the semi- 
amplitude of the eUipsoidal modulation can be approximated a^ 



(4) 



^Fellipsoidal _ . _ (15 + i/i ) (1 + Ti ) y . 2 . 

=0.15^ -\ — q^mi 

F 3-wi \a J 

(lMorris&Naftilanlll993l : IZucker et alJl2007h in which q = M2 /Mi 
is the mass ratio of the two binary components, ui the primary's 
linear limb darkening coefficient and Ti its gravity darkening coef- 
ficient. 

With Ri/a determined from the eclipse duration (Eq.O, q is 
then known as a function of /. In the case of KOI-74, the primary 
is not in corotation jjut is a rapid rotator, as already suggested by 
Ivan Kerkwiik et al ] (I2OIQI) and confirmed by all measurements of 
the rotational velocity (see Section [2]). Van Kerkwijk et al. (2010) 
argue that rapid rotation has a very significant impact on the el- 
lipsoidal modulation amplitude. Kruszewski ( 1963) presented an 
expression for the Roche lobe potentials that accounts for the ef- 
fects of asynchronous rotation. Although we will use this treatment 
in our light curve models, we will be careful when interpreting the 
mass ratio derived from our light curve modelling effort presented 
in Section [3l4l Our modelling setup is such that the mass ratio is 
the only parameter that is only constrained by the ellipsoidal mod- 
ulation amplitude. Therefore, if the expression for the ellipsoidal 
modulation is invalid (e.g. because it assumes instantaneous ad- 
justment of the star's surface to the ever-changing potential in the 
asynchronous case), this will manifest itself in the mass ratio being 
off. 

The detection of Doppler beaming in the extremely accurate 
Kepler light curve (see van Kerkwijk et al. 2010) gives access to 
a second, independent, way to determine ^ as a function of /. For 
velocities much lower than the speed of light, the Doppler beaming 
amplitude is proportional to the radial velocity amplitude of the A- 
star: 



(Eq. 2 in lBloemen et al.ll201ll) . The beaming factor 



(5) 



5 = 5 + dlnF;,/dln;i (6) 

(iLoeb & Gaudill2003h depends on the wavelength X of the observa- 
tion and on the spectrum of the star. We determine this factor for 
the case of KOI-74 in Section[33l 

With Ml given by the spectral type derived from spectroscopy 
(Section 12]), and Ki derived from the Doppler beaming amplitude 
(Eq. [5]) or spectroscopy (Section |2]), we get the following relation 
where q and / are the only unknowns: 
^3 



— sm I. 



(7) 



This equation allows one to derive the mass ratio without relying 
on the amplitude of the ellipsoidal modulation (as long as the de- 
formation of the primary is small enough not to affect the spectro- 
scopic parameter determinations). The mass ratio q derived from 
the radial velocity amplitude using Eq. [7] is not influenced by the 
rapid rotation, contrary to the mass ratio that can be derived from 
the ellipsoidal modulat i on ani plitude via Eq.|4l 

Recently, iKaplanI feOlOl) showed that one can use an effect 
similar to the R0mer delay to derive the radial velocity amplitude 
of the secondary, in double white dwarf binaries with circular orbits 



^ Van Kerkwijk et al. (2010) use a similar equation but with a sin^ / term, 
which should be sin^ i. 



and mass ratios significantly different from 1. In such systems, a 
light travel time difference causes the time between the primary 
and secondary eclipses to be different from P/2 by 

bt = (K2-Ki)— (8) 
nc 

in which c denotes the speed of light (see lKaplanll20 1 Ol) . After sub- 
stitution of ^ = fi^^ following expression: 



ncAt\' 



PKi 



■ 



(9) 



If we can measure this R0mer delay, it would allow us to derive 
the mass ratio in yet another way, this time independently from the 
mass of the primary which was estimated from spectral analysis. As 
far as we are aware, this technique has not yet been applied to any 
system in practice. KOI-74 is a good candidate to put the theory 
to the test, since the expected time difference is of the order of a 
minute. 



3.3 Gravity darkening, limb darkening and Doppler 
beaming coefficients 

The gravity darkening coefficient (see e.g. IClaretl l2003h of the 
primary was calculated by integrating ATLAS model spectra 
JCastelli & Kuru cz 2004) over the Kepler bandpass. We took into 
account the estimated reddening of E{B — V) = 0.\5 (Kepler Input 
Catalo g) by reddening the model spectra following ICardelli et akl 
Reddening marginally influences the gravity darkening co- 
efficient because it is bandpass dependent. Assuming solar metal- 
licity, vturb = 2 km/s, 7;ff = 9500zb250K an d log(g) = 4.3 ± 0.1, 
and using Eq. 1 from iBloemen et al.1 (l201lh . we found the grav- 
ity darkening coefficient to be = 0.55 ± 0.05. We have used a 
physical gravity darkening coefficient of dlogT/dlogg = 0.25. 

Using the same assumptions, we computed limb darkening co- 
efficients for th e A- sta r . We a dopted the 4-parameter limb darken- 
ing relation of IClareJ (|2004 equation 5) with ai = 0.576, ^2 = 
0.118, as = -0.039 and (24 = -0.016. For the white dwarf com- 
panion, we used a model atmosphere for a DA white dwarf with 
Teff = 13000K and \ogg = 6.5 ( Koeste3 |2oTqI) and found ai = 
0.372, a2 = 0.518, = -0.540 and (24 = 0.178. 

For the broadband Kepler photometry, the wavelength- specific 
beaming factor (B in Eq.|5]) has to be replaced by a bandpass- 
integrated photon weighted beaming factor 

Je^XF^BdX 



(B)- 



(10) 



(Eq. 3 in iBloemen et al.ll201lh in which is the response func- 
tion of the Kepler bandpass and B the monochromatic beaming 
factor (Eq. [6]). Taking reddening into account, the beaming fac- 
tor is found to be {B) = 2.19 ± 0.04, which compares well with 
Ivan Kerkwiik et allbOlOh 's value of 2.21. For an unreddened spec- 
trum, the beaming factor would be (B) = 2.22± 0.04. 



3.4 Modelling code and MCMC setup 

Our light curve modelling code, LCURVE (for a description of the 
code, see Copperwheat et al. 2010, Appendix A), uses grids of 
points on the two binary components to calculate the total flux that 
is visible from the system at different orbital phases. It accounts 
for Doppler beaming, eclipses and transits, ellipsoidal modulation 
and reflection effects. The code can also account for lensing effects 
(implemented following iMarsh>2001.) , which occur when the white 
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Figure 4. Illustration of the grids of points we used on the two binary com- 
ponents to model the light curve of KOI-74 (only about one point in 4 is 
shown on the A-star, and about one in 10 on the white dwarf). To achieve a 
high enough numerical precision, at the transit phases, a denser strip is used 
on the primary at the location where it gets occulted by the white dwarf. 
The system is shown with an inclination of 86.4°, at orbital phase 0.02. 

dwarf transits the primary. In the models of KOI-74 presented here, 
however, we did not include these lensing effects. The white dwarf 
radius we would find by including lensing would be slightly larger 
but we estimate the difference at about one percent only, which is 
far lower than the uncertainty on the parameter. In addition, esti- 
mation of the lensing requires knowledge of the mass of the white 
dwarf, which owing to the difficulty in deducing the mass ratio 
from the ellipsoidal variations, was not easily calculated during the 
MCMC runs. 

A typical light curve fit is shown in Fig. [3] together with the 
phase folded light curve and the residuals. Because of the large 
difference in the radii of the stars, it turned out to be difficult to 
get the numerical noise at a lower level than the scatter on the ob- 
servational datapoints of the superb Kepler light curve. The grid 
on the A-star would have to consist of millions of points to get to 
the required accuracy level, which is particularly difficult during 
the transits. We therefore replaced the grid of uniformly distributed 
points, at the transit phases, by a grid with a denser strip on the 
A-star at the region of the star that gets occulted during the white 
dwarf transits. A graphical representation of the grid is shown in 
Fig.m We used 92 716 points (on 270 latitude strips) on the A-star 
outside the transits and 178 468 points during the transits (adding 
7 latitude and 4 longitude points per coarse grid point in the strip). 
The flux from the white dwarf was modelled with a 12 724-points 
grid (100 latitude strips) at all orbital phases. 

To account for a possible detection of the R0mer delay, we 
introduced a 5^ parameter and computed the light curve at phase ^' 
instead of the requested phase (|), with 

<^'^<^+|(^-0.5). (11) 

This essentially results in a maximum R0mer delay of —ht at the 
phase of the eclipse, while the phase of the transit is untouched. 



In an attempt to correctly model the ellipsoidal modulation 
taking into account the effects of rapid rotati on, we imp lemented 
the adapted Roche lobe potentials as given bv lKruszewski (1963J. 
The modelling procedure, using Markov C hain Monte Car l o sim- 
ulations, is identical to the one used in iBloemen et al 
Sections 3.3 and 3.4) for the modelling of the light curve of 
KPD 1946-^4340. In the initial runs using the full dataset, the or- 
bital period (Porb)' the time zero point (ro), the inclination (/), the 
temperature of the white dwarf (72), the mass ratio (^), the beam- 
ing factor {{B)) and the R0mer delay (50 were kept as free param- 
eters. The computation of one synthetic light curve at the 28 086 
time points of the dataset, for which we computed the light curve at 
about 127 000 phases to be able to account for the finite integration 
time of the observations, took about 5 minutes of CPU time. Due to 
the strong degeneracy between the inclination and the mass ratio, 
the MCMC chains took too much time to sample the whole param- 
eter space. We therefore fixed the orbital period (Porb) ^nd the time 
zero point (ro) at the optimal values of the initial runs, and used 
the phasefolded light curve (see Section lTTI for details) for the final 
runs. 

To account for the finite integration times of the observations, 
we oversampled our light curves in time space. The 30 s phasebins 
during eclipses and transits were oversampled by a factor 5, and the 
1800 s phasebins out of eclipse by a factor 3. For the oversampling, 
we approximated the effective integration time for each bin by 

(7)^ + (x)^, in which 7 is the integration time of the binned short 
cadence data points (about 58 s) and x is the width of the bin (30 s 
or 1800 s). To explore the effects of rapid rotation, we ran MCMC 
chains assuming corotation (as in the case of KPD 1946-1-4340), 
as well as chains treating the spin rate as a free parameter with 
vsin/ = 150 ± 10km s~^ as a prior. This prior was implemented as 
a constraint on (27i/?/P) sin/. In total, we computed over 450000 
light curves in the corotation chain, plus about 400 000 light curves 
in the fast rotation chain with spectroscopic prior on the spin rate. 
Every tenth model that was computed was stored and used to deter- 
mine the system parameters. 

We set the limb darkening coefficients to the values found 
in Section 13.31 and used priors on the flux weighted tempera- 
ture of the primary {Ti = 9500 ± 250 K, taken from spectroscopy, 
see Section |2]), the gravity darkening coefficient of the primary 
(Pk = 0.55 ± 0.05, see Section 13.31 ) and the radial velocity am- 
plitude of the primary (Ki = 15.4 ± 0.3 km s~^ taken from spec- 
troscopy, see Section [2]). If the amplitude of the Doppler beaming 
effect is consistent with t he radial velocity anip litude, as was the 
case for KPD 1946-^4340 (iBloemen et aDl201lh . the results of the 
analysis will be identical when the prior on Ki would be replaced 
by a prior on (B) . In our approach, putting a prior on Ki and treating 
(B) as a free parameter, we can judge whether the Doppler beam- 
ing effect has the expected amplitude by comparing the (B) that 
results from our MCMC analysis with the one found from atmos- 
phere models in Section [331 



3.5 MCMC results and discussion 

As can be seen in Fig. [3] the light curve model fits the observed data 
very well. There is some structure in the residuals around the transit 
ingress and egress (around orbital phase 0.5), but it is barely signif- 
icant. This structure can for example result from a small difference 
between the assumed limb darkening coefficients (which were de- 
rived for a spherical star) and the true limb darkening of the A-star. 
The parameters derived from our MCMC runs that account for the 
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Table 2. Properties of KOI-74. The mass of the primary, the effective tem- 
perature and the rotational velocity are derived from spectroscopy. The 
other parameters are the result of our MCMC analysis of the Kepler data. 
The two different mass ratios are derived from the information contained 
in the ellipsoidal modulation amplitude (^eii) and the radial velocity (^rv). 
The mass of the secondary is derived using ^rv. 

Primary (A- star) Secondary (WD) 
/'orb (d) 5.188675(4) 
/(deg) 87.0 ±0.4 

Teff (K) 9500 ± 250 14 500 ± 500 

^ell,corot 0.052 ±0.004 

^ell,vsini 0.047 ±0.004 

^RV 0.104 ±0.004 

R{^q) 2.14±0.08 0.044±0.002 

M (M©) 2.2 ± 0.2 0.228 ± 0.014 

vsin/(kms-i) 150±10 

5fR0merW -56 ±17 

{B) 2.24 ±0.05 



, using the prior on v sin /, are summa- 
rized in Tablel2l lRowe et al.l(l201Ql) determined Mi = 2.2 ± O.2M0 
based on the spectra l type. Using the evolutionary model grids of 
iBriquet et alJ (l201lh and our spectroscopic determinations of T^ff 
and logg, we find the same result, which we adopted for our anal- 
ysis. The results of the MCMC chains for the corotating case are 
nearly identical except for the derived mass ratios. The mass ratio 
from the models assuming corotation is also listed in the table, to 
allow for comparison with previously published values. The most 
important difference between our results and those previously pub- 
lished, is that our inc lination / = 87.0 ± 0.4° differs by ~ 3a from 
the value derived by iRowe et all (l2010h . / = 88.8 ± 0.5°. In that 
paper, the inclination was determined by fitting th e eclipses and 
transits using the analytical formulae of Mandel & Agol (20o3), 
which include limb darkening but not gravity darkening. Further- 
more, it is possible that the authors did not account for the finite 
exposure times, which significantly smear out the eclipse ingresses 
and egresses in the long cadence data t hey had avai l able. Fitting 
only the 43 d of long cadence data that iRowe et aP (l2010h used, 
we find an uncertainty on the inclination of 0.9°, compared to their 
more optimistic value of 0.5°. The 4 months of short cadence ob- 
servations allowed us to determine the inclination more reliably. 

The distribution of the inclination values in our MCMC runs 
using the prior on vsin/ is shown in the bottom panel of Fig. [5] 
The two top panels of the figure show our radius estimates (relative 
to the separation of the two stars, a) as a function of the inclina- 
tion. The blue points with error bars indicate the values found by 
Ivan Kerkwiik et al who did not derive the inclination inde- 

pendently but instead adopted the value of iRowe et al ] (l2010h .Due 
to our lower preferred inclination, the radii in ferred in our analy- 
sis are slightly larger than the ones obtained bv lvan Kerkwiik et al] 
(l2010h . 

The third panel of Fig. |5] shows the mass ratio as a function 
of the inclination. The red line is the theoretical relation, assuming 
corotation, based on the constraints offered by the transit duration 
(Eq.[3]) and the ellipsoidal modulation amplitude (Eq.|4]). The dis- 
tribution of points from our corotation runs, falls nicely around this 
line (not shown). The blue line shows the mass ratio derived from 
the spectroscopic Ki , or equivalently, as we will see further in this 
Section, the Doppler beaming amplitude. The discrepancy between 
the mass ratios derived from the ellipsoidal modulation amplitude 
on one hand, and from the radial velocity information on the other 




86 87 88 89 
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Figure 5. Illustration of the correlation between the inclination of the sys- 
tem and the mass ratio and radii of the stars. The bottom panel shows the 
distribution of the models that were accepted in our MCMC runs that ac- 
count for the rapid rotation of the primary. The dashed (dot- dashed) line 
indicates the 68 (95) per cent confidence interval. The black dots represent 
a random selection of the models. The blue dots with error bars show the 
results from van Kerkwiik et al. (2010). The red line on the mass ratio plot 
is the theoretical relation based on Eq. [3] and Eq. |4] assuming corotation. 
The blue line shows the mass ratio derived from K\ (Eq.[7j. 



hand, was already noted by[ van Kerkwijk et al The results 

of the models that account for rapid rotation of the primary are plot- 
ted with black dots. It is striking that these models result in lower 
rather than higher mass ratios, thus only increasing the discrepancy 
compared to the models assuming corotation. The lower preferred 
inclination value also contributes to the increase in the discrepancy 
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Figure 6. Distribution of the R0mer delay as fitted by our MCMC runs. The 
delays expected for mass ratios of g = 0.050 and q = 0.104 are indicated 
for comparison. 
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Figure 7. Time-frequency analysis of the residuals of the short cadence data 
after subtracting our best fitting binary light curve model. The observed 
variability at ~ 1.3 d~^ might be related to the spin period of the primary 
star. 



between the two mass ratios, compared to the discrepancy shown 
by[ van Kerkwiik et all (l20ld) . This can easily be understood: as / 
gets lower, R\/a has to increase to make the model fit the observed 
transit duration (Eq. [3]); but a higher Ri/a would lead to a higher 
ellipsoidal modulation amplitude, which is then compensated in the 
simulations by lowering q^w (Eq.|4l). 

The beaming factor derived from the MCMC chains, {B) = 
2.24 ± 0.05, agrees with the beaming factor derived from spec- 
troscopy in Section [33l f(^) = 2. 19 ±0.04, taking reddening into 
account). We could thus have modelled the binary equally well 
without spectroscopic radial velocity information, just relying on 
the Doppler beaming amp litude in the Kepler light curve, verifying 
Ivan Kerkwiik et alj (|201Q|) 's approach. 

One can argue that the mass ratios determined from the el- 
lipsoidal modulation amplitude and radial velocity information can 
be brought into agreement by assuming a high contamination of the 
Kepler light curve by background stars, since this would increase 
the observed ellipsoidal modulation amplitude. The consistency be- 
tween the observed Doppler beaming amplitude (which would also 
increase if one assumes a higher contamination) and the spectro- 
scopic radial velocity amplitude, however, implies that our assump- 
tion of no background contamination in the Kepler light curve has 
to be correct up to a few per cent. 

We were also able to measure the R0mer delay at 5^ = — 56 ± 
17 s, as can be seen on the distribution plot from our MCMC 
chains on Fig. [6] The expected R0mer delay for a mass ratio of 
q = 0.104 is bt = -63.2± 1.2 s, while for q = 0.050 we would ex- 
pect 5? = — 139 ± 3 s (the error bars account for the uncertainty on 
Ki measured from spectroscopy). The R0mer delay estimates only 
depend on Poj-b^ the mass ratio. With Poj-^ ^nd Ki firmly 

established from both spectroscopy and photometry, the measure- 
ment of the delay proves that the true mass ratio of the system is 
the one derived from the radial velocity amplitude, under the con- 
dition that the orbit is circular. Note that, while Ki can be derived 
from the Doppler beaming amplitude, a measurement of the R0mer 
delay allows one to also derive K2 directly from the light curve if 
Ml is known. 



3.6 Variability in residuals 

We performed a time-frequency analysis on the residuals of the 
short cadence data after subtracting our best model. The time- 
frequency diagram is shown in Fig. [71 We detected sign ificant vari- 
ability (following the criteria given in lDegroote et al. 12009), but all 
at low amplitudes (~ 10 parts per million). We found variability 
with periods of ~ 3 d, which is a known instrumental artefact (see 



[Christiansen et al.l201lh . We also found variability with a period of 
0.5918 ±0.0015 d, with an amplitude that changes in time (see Fig. 
IT}. Van Kerkwijk et al. (2010) suggested that this signal could be as- 
sociated with the spin period of the A- star. Given our spectroscopic 
value for vsin/ and the radius determination from our light curve 
analysis, we expect the spin period of the A-star to be 0.72 ± 0.06 d, 
which differs by 2a from the detected periodicity. 

We detected additional variability of similar amplitude in the 
long cadence Q2 data at frequencies between 15 and 20 d~^. We 
could find a similar signal in the observations of KIC 6889190, 
which is observed on the same CCD close to KOI-74, and it can be 
removed su ccesfully using the cotrending basis vectors technique 
described in [Christiansen et aP (l201lh (Tom Barclay, priv. comm.). 
This confirms that it is instrumental in nature. 

We did not find significant residual power at the orbital period 
which proves that the LCURVE model fits the data very well. 



4 SUMMARY 

We have analysed 212 d of Kepler data (QO, Ql, Q2 and Q3) of 
the eclipsing binary KOI-74, as well as supporting spectroscopic 
observations. We modelled the light curve using the LCURVE code, 
accounting for ellipsoidal modulation, reflection effects, Doppler 
beaming, R0mer delay, eclipses and transits. Using Markov Chain 
Monte Carlo simulations, we determined various system param- 
eters of KOI-74, which are summarised in Table [2] We find a 
lower orbital inclination of / = 87.0 ± 0.4 ° compared to the dis- 
covery pape r and first analysis p r esente d bv lRowe et al ] (l2010h and 
adopted by Ivan Kerkwiik et alJ tOld) , i = 88.8 ±0.5°. The dif- 
ference propagates to our values of other parameters such as the 
radii. It lowers the mass ratio derived from the ellipsoidal mod- 
ulation amplitude to q = 0.052 ± 0.004 assuming corotation or 
q = 0.047 ± 0.004 accounting for a rapidly rotating primary (us- 
ing vsin/ = 150zblOkms~^ as a prior). 

The amplitude of the observed Doppler beaming, which 
makes the primary become brighter when the star moves in the 
direction of the Kepler satellite in its orbit, is in perfect agree- 
ment with what is expected from the spectral type and radial ve- 
locity amplitude of the star, which we determined to be Ki = 
15.4 ± 0.3 km s~^ from spectroscopy. From the primary's mass of 
2.2±O.2M0 the mass ratio derived from the radial velocity is 
^ = 0.104 ±0.004. 

We also report the first detection of R0mer delay in a light 
curve of a compact binary. This delay, which amounts to 56 ± 17 s, 
is exactly as long as one would expect for a mass ratio of ^ ~ 0.1, 
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and is in contradiction with the lower mass ratio derived from the 
eUipsoidal modulation amplitude. Van Kerkwijk et al. (2010) pre- 
ferred the mass ratio derived from the Doppler beaming amplitude 
over the mass ratio derived from the ellipsoidal modulation am- 
plitude. We find that the spectroscopic radial velocity amplitude 
as well as the R0mer delay leave no doubt that this was indeed the 
correct assumption, and that the secondary of KOI-74 is a low mass 
white dwarf. 

As a result of our lower preferred orbital inclination value, the 
discrepancy between the mass ratio determined from the ellipsoidal 
modulation amplitude and the higher mass ratio determined from 
ra dial velocity information or R0mer delay has increased compared 
to Ivan Kerkwijk et aP bOld) , and now amounts to a factor 2. Our 
attempt to account for the effect of the rapid rotation of the primary, 
increased the discrepancy even further. Our results imply that one 
has to be very cautious when adopting mass ratio estimates derived 
from the ellipsoidal modulation amplitude, in particular if there is 
no firm proof of corotation. We are not aware of any theoretical 
explanation of the reduction in ellipsoidal modulation as a result of 
rapid asynchronous rotation. 
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